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A linear filter to reconstruct the ISW effect from 
CMB and LSS observations 
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Abstract — The extraction of a signal from some observational 
data sets that contain different contaminant emissions, often at 
a greater level than the signal itself, is a common problem in 
Astrophysics and Cosmology. The signal can be recovered, for 
instance, using a simple Wiener filter. However, in certain cases, 
additional information may also be available, such as a second 
observation which correlates to a certain level with the sought 
signal. In order to improve the quality of the reconstruction, it 
would be useful to include as well this additional information. 
Under these circumstances, we have constructed a linear filter, 
the linear covariance-based filter, that extracts the signal from 
the data but takes also into account the correlation with the 
second observation. To illustrate the performance of the method, 
we present a simple application to reconstruct the so-called 
Integrated Sachs-Wolfe effect from simulated observations of the 
Cosmic Microwave Background and of catalogues of galaxies. 



I. Introduction 

THE issues of signal reconstruction and component sep- 
aration have become of major importance in the field 
of Astrophysics and Cosmology, covering a large range of 
applications. In particular, Cosmic Microwave Background 
(CMB) observations require the development of sophisticated 
image processing techniques, in most cases on the sphere, to 
extract all the valuable information encoded in this type of 
data. The CMB is a relic radiation that was emitted shortly 
after the Big Bang, when the universe was ~ 380000 years 
old and has travelled since then - almost unhindered - towards 
us, carrying a wealth of information about the origin and 
evolution of the universe (see e.g. [1], [2] and references 
therein). From the standard theory of inflation, the statistical 
distribution of the CMB anisotropies is expected to follow an 
isotropic and Gaussian random field on the sphere. Different 
physical phenomena leave their imprint in the CMB radiation 
that we observe today. In particular, a relevant contribution 
is given by the so-called Integrated Sachs-Wolfe (ISW) effect 
[3], which is due to the evolution in time of the gravitational 
potential, in the linear regime, of the Large Scale Structure 
(LSS) and its interaction with the CMB. This effect is of 
great interest since it provides a direct indication of either 
the presence of dark energy in the case of a flat universe or 
the existence of spatial curvature [4]. 

Using only CMB observations, it is difficult to obtain a 
direct detection of the ISW signal, since its contribution to 
the total CMB signal is, in general, small, and is difficult to 
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disentangle it from other physical effects that also constitute 
the observed CMB picture. However, given the origin of this 
effect, a spatial correlation between the ISW and the LSS is 
expected. Therefore, it was suggested [5] that the signal of 
the ISW effect could be statistically detected by looking for 
cross-correlations between CMB and LSS observations. Fol- 
lowing these ideas, different groups have obtained a statistical 
detection of the ISW effect using the CMB data provided by 
the WMAP NASA satellite and different LSS surveys at levels 
between 2.5er and 4.5cr [6]— [16]. 

Nonetheless, it would be desirable to obtain not only a 
statistical detection of the ISW signal but also the actual 
map of this effect. A possible procedure would be to apply 
a reconstruction image technique (e.g. Wiener filter [17], 
maximum-entropy method [18]) to the CMB map in order 
to extract the ISW effect from the rest of contributions of 
the primordial radiation. However, given the weak level of 
the signal such a simple procedure is not expected to provide 
good results. Also, techniques that make use of multifrequency 
observations (see [19] and references therein) are not useful 
since the signal has the same frequency dependence as the 
other dominant CMB contributions. Thus, instead, we propose 
to recover the ISW map applying a filter, that we have called 
linear covariance-based (LCB) filter, to the CMB data that 
takes also into account the information coming from the cross- 
correlation between the CMB and the LSS surveys. 

The outline of the paper is as follows. Section [II] presents 
our method to reconstruct a signal embedded in a noisy back- 
ground and correlated with a second observation. Section [TlTl 
shows the performance of the proposed LCB filter in a simple 
application, the recovery of the ISW signal from simulated 
data sets. The robustness of the methodology against certain 
non-idealities is studied in section [TV] Finally, we present our 
conclusions in section [V] 

II. The method 

A random signal a measured on the sphere, such as the 
CMB, is usually expanded in terms of the spherical harmonics 
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where 6 £ [0, %] and <f> € [0,27r] stand for the colatitude and 
longitude in spherical coordinates, respectively, and ai m are 
the coefficients of the expansion. For a statistically isotropic 
signal, the variance of the coefficients is independent of m: 
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where the averages are to be taken over statistical ensembles 
and the CV's constitute the angular power spectrum. For the 
case of an isotropic Gaussian random field on the sphere, the 
power spectrum completely characterises the signal. 

A. Wiener filter (WF) 

Let us assume that we have an observation d that contains a 
certain signal of interest s plus some generic noise n, both of 
which are assumed to be Gaussian and isotropic. One possible 
way to obtain an estimation s of the signal is to apply the 
classical WF to the data, which, at each (£, m) mode, is given 
by: 

C s 

* em = c s + c n dtm (3) 

where C| and Cf correspond to the power spectrum of signal 
and noise respectively. It can be shown that WF is the linear 
filter that minimises the variance of the reconstruction error. 

It is well known that the power spectrum of the WF 
reconstruction is biased towards values lower than the true 
signal, with the bias depending on the signal-to-noise ratio 
of the data. In particular, from the previous equation, it is 
straightforward to show that the expected value of the power 
spectrum for the reconstructed signal is given by 

(Off 



cs + cv 



(4) 



or, alternatively, we can write the expected bias of the signal 

as 

(b) = CI (Ci) = (5) 

Note that in the limit of infinite signal-to-noise ratio, the WF 
goes to unity and the bias becomes negligible. 

B. Linear covariance-based (LCB) filter 

Let us assume that we have an additional second observation 
g which is also isotropic and Gaussian and that is correlated 
with the signal s. It would be desirable to include the informa- 
tion of g in order to reconstruct s. The signals are characterised 
by their angular power spectra and also by their cross power 
spectrum: 

Cf = (ge m s* e , m ,) (6) 

which we assume to be (statistically) knowrQ. The covariance 
matrix C(£) of g and s at each multipole £ is constituted by the 
expected power spectra of the signals (diagonal elements) and 
their cross-power spectrum (off-diagonal elements). It is useful 
to calculate the Cholesky decomposition of the covariance 
matrix, which satisfies C(£) = L(£)L T (£), where L is a lower 
triangular matrix. It can be shown that gi, n and si m can be 
written as 



Ln(£) 
L l2 {£) L 22 {£) 
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(7) 



1 In the particular case of the IS W recovery, this corresponds to the expected 
cross-correlation between the ISW signal and the projected galaxy density. 
This quantity is given in terms of cosmological parameters that define the 
expansion and the energy content of the Universe, and the astrophysical laws 
that guide the LSS evolution. 



where h and j are uncorrected Gaussian variables of unit 
variance, sometimes referred to as hidden variables. Note that 
the elements of L{£) relate to the elements of the covari- 
ance matrix C(£) as L\\ = \JCf,L\i — Cf s / \[Cf and 



J i 

L22 = y/\C(£)\ JCj, where \C(£)\ is the determinant of the 
covariance matrix at each £ mode. 

From the previous equation, we can trivially write s as a 
function of g and j as 

Li 2 (£) 
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Let us recall the observational data d, which can now be 
written as 



dim — Sim + nim — 
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In order to recover s, we can now proceed in the following 
way. The first term of equation ([8j can be easily calculated 
from g and the covariance matrix of the signals. The second 
term can be estimated by applying WF to the modified data d 
given by 
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Therefore, the final estimation s, at each multipole, is obtained 

as 
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This simple procedure allows one to take into account the 
second correlated observation g in the reconstruction of s. 
This improves the quality of the reconstructed signal that also 
has the right correlation with g. 

It can be easily shown that the expected value of the power 
spectrum for the reconstructed signal is given by 

/rS s_ {Cr) 2 {\C{£)\+CjCf) + \C{£)\ 2 

Ki) C 9 e (\C(£)\ + C 9 e C?) ( ' 

As we are using WF to recover a part of the signal, the power 
spectrum of the reconstruction is again biased. However, it can 
be shown that this bias is always lower or equal than the one of 
the WF reconstruction. To prove this statement, we first note 
that the cross power spectrum must satisfy (C^ s ) 2 = aCfCf, 
with < a < 1 . Taking this property into account and making 
use of equation (fl2l i. we find that the bias introduced by the 
LCB filter is given by 



(b)= 



Cf 



l-a 



(13) 



As one would expect, the bias is a function of Cf, Cf and a, 
but does not depend on Cf. It is interesting to point out that 
for a = 0, the LCB filter defaults to WF, since in this case g 
does not contain information about s. For a = 1, the bias goes 
to zero, and the reconstructed signal is simply proportional to 
the observation g. Also, by comparing equations (f5]l and ( fT3l , 
and since < a < 1, it becomes apparent that the bias of the 
LCB filter reconstruction is always lower or equal than the 
one of WF. 
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Finally, we would like to remark that the described proce- 
dure can be easily generalised to include additional correlated 
signals, providing the expected value of the full covariance 
matrix is known. 

III. Reconstructing the ISW map 

In order to test the performance of our method, we have 
considered the problem of reconstructing the ISW map from 
CMB and LSS observations. In our case, the signal s cor- 
responds to the temperature fluctuations caused by the ISW 
effect, the observed data d is a CMB observation, the noise 
n is all the signal present in the CMB map that is not ISW 
and g is a catalogue of galaxies (i.e. a map of the galaxy 
density projected along the line of sight) which is correlated 
with the ISW pattern through the gravitational potential which 
affects to both. The reconstruction of the ISW map is a 
very challenging problem since the signal is very weak, in 
particular the signal-to-noise ratio for the considered example 
(maps with ~ 2-degree resolution) is around 0.4. Although 
real CMB and LSS data are affected by non-idealities, like 
incomplete sky coverage and residual signals coming from 
astrophysical contaminants, in this section we will consider 
full-sky observations free of residual contaminants. The reason 
is that this simpler scenario allows one to understand better the 
performance of the LCB filter. However, in section [IV] we will 
study the robustness of the method in more realistic situations. 

We have generated correlated simulations of CMB and 
catalogues of galaxies making use of equation (0 (see also 
[20]). More specifically, for a given model and at each 
multipole, first we calculate the covariance matrix C(£) (using 
a modified version of CMB FAST [21]) and its corresponding 
Cholesky decomposition L(£). Then, we generate two uncor- 
rected Gaussian variables, hg m and ji m , of zero mean and 
unit dispersion. These quantities are then combined following 
equation ((T) to obtain the coherent ISW and galaxy catalogue 
simulations. Finally, an independent Gaussian CMB simula- 
tion is generated (without including the ISW contribution) and 
added to the ISW map to obtain the simulated data d. Since the 
ISW effect is important at large angular scales, our simulated 
maps have a resolution of ~ 2 degrees, corresponding to the 
HEALPixQ parameter n s id e = 32. For the CMB power spectra 
the flat A-CDM model that best fits the 3-year WMAP data 
[22] has been used. According to this model, the Universe 
has undergone a recent epoch of accelerated expansion, and 
the energy causing this acceleration amounts, at the present 
time, around 75 per cent of the total energy budget of the 
Universe. We have neither included a beam nor instrumental 
noise, since, for the current and future CMB observations, then- 
effect is negligible at the considered angular scales. For the 
galaxy catalogue, we have generated several idealised surveys 
centred at different redshifts z (an astrophysical parameter 
related to the distance of the galaxies from the observer). It is 
known that galaxies are biased tracers of the underlying matter 
density field present in the Universe, in the sense that they 
tend to populate most of the times the most overdense regions 

2 HEALPix is a hierarchical equal area isolatitude pixelization of the sphere. 
For more information see http://healpix.jpl.nasa.gov/ 



of the Universe. This bias is dependent of the galaxy type, 
and a priori may depend as well on the age of the Universe 
(also related to the redshift parameter z). For simplicity we 
have considered a constant bias and a population of galaxies 
with no intrinsic evolution with redshift. It can be shown 
from equation (TTTb that the reconstruction of the signal is 
independent of a constant bias, since it acts as a global 
normalisation of the galaxy catalogue. The cross-power spectra 
between the CMB and the different catalogues of galaxies have 
been obtained following linear theory within the frame of the 
concordance ACDM cosmological model. This assumes that 
perturbations in the density field and gravitational potential 
are small, which is appropriate for the large scales relevant 
to the ISW effect (we know that the universe is isotropic 
and homogeneous in the very large scales). This computation 
traces the freeze-out of the growth of the gravitational potential 
wells due to the onset of an accelerated expansion in the 
Universe, and couples it to the gravitational energy that CMB 
photons gain or lose when crossing those potential wells. 

As a measurement of the quality of the recovered signal, 
we will use the correlation coefficient p between the input 
and reconstructed maps which is defined as 



cr s a$ 

where a s and a§ are the dispersion of the input and recon- 
structed map respectively. Values of p closer to one provide 
better reconstructions. As an illustration, Fig. 03 shows a 
simulation of the ISW, the CMB (including the ISW) and a 
catalogue of galaxies centred at a comoving radial distance r 
corresponding to a redshift z = 0.8. The width of the radial 
shell containing the galaxies is roughly O.lr. Fig. [2] shows 
the reconstructed ISW map using the LCB filter (top panel) 
for the set of simulations given in Fig. [T] For comparison, 
the WF reconstruction is also shown (bottom panel). It is 
apparent that the map recovered with the LCB filter follows 
better the general structure of the input ISW map compared to 
that obtained with WF. This can also be seen from the cross- 
correlation coefficient between the input and reconstructed 
map which is 0.88 and 0.63 for the LCB and WF, respectively. 

For a given CMB observation, the quality of the recovered 
ISW map obtained with the LCB filter could be significantly 
different depending on the catalogue of galaxies chosen for the 
reconstruction, since the cross-correlation between the ISW 
and the LSS is redshift dependent (as well as depends on the 
particular type of galaxies, the completeness of the survey, 
etc.). To provide an insight of the importance of this effect on 
the quality of the reconstruction, Fig. [3] shows the correlation 
coefficient (averaged over 1000 simulations) for the LCB filter 
using a CMB map and different LSS surveys, centred at 
redshift z (black solid thick line). We have also considered two 
widths in the range of redshifts observed by the survey: a wider 
bin centred on z (or on distance r) with a width of O.lr (left 
panel) and a narrower bin of width 0.05r (right panel). The 
thinner solid curves provide the la region around the mean 
p obtained from the simulations. For comparison, the results 
for the reconstruction from WF (red dashed thick line) and 
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Fig. 3. Mean correlation coefficient for the ISW reconstruction using the LCB filter obtained from 1000 simulations (black solid thick line) and the 
corresponding lcr region. Two cases are considered: a wider bin in z (left panel) and a narrower bin (right panel). For comparison the mean correlation 
coefficient and the corresponding lcr region for the WF reconstruction are also shown (red dashed lines). 




Fig. 1. The figure shows coherent simulations of an ISW map (top 
panel), a CMB map including the ISW effect (middle panel) and a catalogue 
of galaxies (bottom panel). The catalogue of galaxies has arbitrary units 
(the reconstruction is independent of the normalisation) and is centred at a 
comoving radial distance r corresponding to a redshift z = 0.8. The width 
of the radial shell containing the galaxies is approximately O.lr. The units of 
the CMB and ISW maps are microkelvin. 




Fig. 2. Reconstruction of the ISW signal from the simulated data of 
Fig. [2 using the LCB filter (top) and the WF (bottom). The figures should 
be compared to the input ISW (top panel of Fig. [TJ. The cross-correlation 
coefficient p is equal to 0.88 and 0.63 for the LCB and WF reconstructions, 
respectively. 

the corresponding lcr region (red dashed thin lines) are also 
shown. Note that the WF reconstruction is independent of the 
considered LSS survey, since it does not use this information 
(the small deviations from a constant line are due to the finite 
number of simulations). From Fig.[3]it is apparent that the LCB 
filter outperforms the WF in all the considered cases. One can 
also notice that the surveys with a wider bin in z work better 
than those with a narrower bin. In particular, for the considered 
cases, we find that the best reconstruction is obtained for a 
survey centred in z = 0.8 for the case of the wider bin. This 
indicates that having a larger number of observed galaxies 



JOURNAL OF MeX CLASS FILES - 



5 



10 3 




10 1 1 1 1 1 1 ' 

5 10 15 20 25 30 

Multipole (I) 



Fig. 4. Mean power spectrum of the ISW reconstruction using the LCB filter 
obtained from 1000 simulations (black squares) compared to the input model 
(blue solid line). The blue dotted lines indicate the Icr region allowed by the 
cosmic variance and the error bars give the dispersion of the reconstructed 
power spectrum obtained from simulations. The theoretical expectation for 
the power spectrum of the reconstruction is also plotted (black dashed line). 
For comparison, the same quantities are given for the Wiener filter case (red 
circles and red dashed line). The reconstructions correspond to the wider bin 
case, centred at z = 0.8. 

increases the correlation between the ISW and the considered 
survey and that, therefore, the reconstruction of the ISW signal 
is improved. However, we should remark that the bias between 
the galaxy and the underlying matter distribution is more 
unlikely to remain constant (as we have assumed) for wide 
bins. For this reason, a more detailed study of the bias within 
the redshift/distance bin would be required. 

It is also interesting to point out that the quality of the 
reconstruction is higher in the hottest and coldest spots of the 
ISW map, since the signal is comparatively larger in those 
regions. For instance, for the best reconstruction case, the 
average correlation coefficient increases from p = 0.85, when 
all the pixels are considered, to p = 0.94, when only the pixels 
outside the la region are considered. 

We have also studied the power spectrum of the recovered 
ISW map, which provides an alternative measurement of the 
quality of the reconstruction. As mentioned in section [TT1 we 
expect the power spectrum of our recovered ISW map to be 
biased at a certain level. Fig. [4] shows the power spectrum 
of the reconstructed ISW using the LCB filter averaged over 
1000 simulations (black squares) for the case of a catalogue 
of galaxies covering a wide bin centred at redshift z = 0.8, 
compared to the input ISW model (blue solid line). The blue 
dotted lines correspond to the theoretical cosmic variance 
while the error bars give the dispersion obtained from sim- 
ulations. For comparison, the expected reconstructed power 
spectrum given by equation (fl2l i is also plotted (black dashed 
line), showing an excellent agreement with the results obtained 
from simulations. Although, as expected, the reconstructed 
power spectrum is biased, taking into account the weakness 
of the ISW effect, a significant part of the power of the 
signal is present in the reconstructed map. In particular, the 
average reconstructed power spectrum lies within the allowed 
la cosmic variance region up to I ~ 15. We have also studied, 



for the same case, the average power spectrum and dispersion 
of the reconstructed ISW map using WF (red circles and error 
bars in Fig. @J, finding a much larger bias, in agreement 
with equation (red dashed line). This shows again that 
combining information from the CMB and LSS, as in the LCB 
filter, is very useful to recover the ISW map. 

IV. Robustness of the method 

In the previous section, we have shown the performance of 
the LCB filter to recover the ISW effect from CMB and LSS 
observations in ideal conditions. However, when dealing with 
real data, they will be affected by different non-idealities, that 
should be taken into account in the analysis. Regarding CMB 
data, one of the most important issues is the contamination 
due to different astrophysical emissions, which is particularly 
strong in the Galactic plane. This region is usually discarded 
from the data, by applying a Galactic mask (i.e. replacing with 
zeroes the most contaminated regions). Regarding catalogues 
of galaxies, they do not usually cover the whole sky and, 
therefore, it is necessary to deal with incomplete sky maps. 

In order to test the robustness of the LCB filter on an in- 
complete sky, we have repeated the analysis from the previous 
section after applying a simple mask to the data that covers 
around 25 per cent of the sky (see top panel of Fig. 0. The 
mask has been constructed as the addition of a band centred 
in the equator with latitude |6| < 8° plus a circle centred in 
the southern ecliptic circle with a radius of 40°. The equatorial 
band mimics a simple Galactic mask from a CMB observation, 
discarding ~ 15 per cent of the sky, which is similar, for 
instance, to the area covered by the WMAP Kp2 mask [23]. 
The circle mimics a region not observed by an incomplete 
galaxy survey, as in the case of the NVSS catalogue [24]. 

The top panel of Fig. [5] shows the reconstruction of the 
ISW map using the LCB filter from the simulated data 
of Fig. 03 after applying the considered mask. The cross- 
correlation coefficient p, obtained outside the masked region, 
is 0.81 to be compared to p=0.88 in the ideal case. A more 
detailed study of the effect of the mask on the quality of the 
reconstruction is given in the left panel of Fig. [6] where the 
mean and dispersion of the correlation coefficient for the ISW 
reconstruction is shown for the LCB filter (black solid lines) 
and for the WF (red dashed lines), considering the case of 
the wider bin in z. This should be compared to the results 
obtained for the same case assuming ideal conditions, given 
in the left panel of Fig. [3] We can see that the presence of 
the mask reduces, although only at the level of a few per 
cent, the correlation between the input and reconstructed maps. 
For instance, for the best reconstruction, corresponding to a 
catalogue of galaxies centred at z — 0.8, the mean value of the 
correlation coefficient is reduced from 0.85, in the ideal case, 
to 0.81, when the mask is present. Therefore, the performance 
of the LCB filter is robust against incomplete sky observations. 
Regarding WF, the mean correlation coefficient is very similar 
to the one obtained in the ideal case. 

Another interesting issue is the robustness of the method 
against uncertainties in the knowledge of the covariance matrix 
C(£) or, for the considered case of the ISW reconstruction, 
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Fig. 6. Mean correlation coefficient for the ISW reconstruction using the LCB filter obtained from 1000 simulations (black solid thick line) and the 
corresponding 1<t region (black thin lines) for different galaxy surveys (considering the case of the wider bin in z). The left panel corresponds to reconstructions 
obtained from data where a mask is applied (discarding around one-fourth of the sky) and the right panel shows the results obtained assuming a wrong 
cosmological model. For comparison, the mean correlation coefficient and the corresponding la region for the WF reconstruction are also shown (red dashed 
lines). 
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Fig. 5. Reconstruction of the ISW signal for the case of an incomplete sky 
(top) and when assuming an erroneous cosmological model (bottom). The 
figures should be compared to the input ISW (top panel of Fig. fT} and also 
to the reconstruction in the ideal case (top panel of Fig. |5J. The grey region 
of the top panel indicates the masked region, excluded from the analysis. The 
cross-correlation coefficient p is equal to 0.81 (top reconstruction) and 0.87 
(bottom reconstruction). 

in the knowledge of the cosmological model. In the previous 
examples, we have used the same cosmological model to 
generate the simulations and to reconstruct the ISW map. 
However, in a more realistic case, the assumed cosmological 
parameters may deviate, at a certain level, from their true 
underlying values. 

Different cosmological data sets (such as CMB, LSS galaxy 
surveys, primordial Big-Bang nucleosynthesis, measurements 



of the Hubble constant or supernovae data) are currently 
placing strong constraints on the cosmological parameters, 
giving rise to a consistent picture of the universe, the so- 
called concordance model. Therefore, given the current con- 
straints, large deviations from the present estimation of the 
cosmological parameters are not expected. In order to test the 
sensitivity of the LCB filter to possible errors in the assumed 
cosmological model, we have produced reconstructions of the 
ISW map using a model different from the one present in 
the simulations (which was given by the A-CDM model that 
best fitted the 3-year WMAP data), i.e, we obtain the recon- 
structions using a wrong covariance matrix. This second model 
corresponds to the parameters given by the A-CDM model that 
best fits the 5-year WMAP data, baryon acoustic oscillations 
in the distribution of galaxies and Type I supernovae data (see 
[25] and references therein). The differences between these 
two models correspond to the level of uncertainties that one 
would expect in our present knowledge of the cosmological 
parameters. In particular, the value of the dark energy density 
ft\, one °f the parameters more closely related to the ISW 
signal, differs in a ~ 5 per cent between both models. 

The bottom panel of Fig. shows the reconstruction of 
the ISW map using the LCB filter, assuming the alternative 
cosmological model, from the simulated data of Fig. [T] Al- 
though, in this case, the amplitude of the reconstructed ISW is 
a bit larger (due to the different assumed cosmological model), 
the correlation coefficient (p — 0.87) is still very similar to 
the value found in the ideal case. In addition, the right panel 
of Fig. [6] shows the mean and dispersion of the correlation 
coefficient for the ISW reconstruction in this case, averaged 
over a large number of simulations and considering the case 
of the wider bin in z. Again, the figure should be compared 
to the same case assuming ideal conditions, given in the left 
panel of Fig. [3] From the comparison, it is apparent that the 
assumption of a different cosmological model does not affect 
significantly the mean correlation coefficient between input 
and reconstructed ISW map neither for the LCB filter nor 
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for the WF. This result shows that our methodology is also 
robust against the possible uncertainties in our knowledge of 
the cosmological parameters. 

V. Conclusions 

We have presented a linear filter, the LCB filter, which 
reconstructs a signal from a confusion background taking 
into account additional information from a second observation 
correlated with the sought signal. To study the performance of 
the technique, we have applied it to simulated CMB and LSS 
observations with the aim of recovering the ISW effect. This 
is a very challenging problem, due to the weakness of the 
ISW signal and, as far as we know, this is the first time that 
a technique is presented that attempts to recover an ISW map 
instead of searching for a statistical detection of the signal. 
We have considered different possibilities for the LSS survey, 
in order to illustrate the dependence of the method on the 
characteristics of the chosen catalogue of galaxies. From the 
considered cases, the best reconstruction (mean p — 0.85) 
is obtained for a survey with a wide bin in z and centred 
in z = 0.8. For comparison, we have also recovered the 
ISW signal from CMB simulated maps using a classical WF, 
finding that the LCB filter provides better reconstructions. This 
confirms the idea that combining information from the CMB 
and the LSS survey is useful to recover the ISW map. We 
have also shown that our methodology is robust when dealing 
with incomplete sky observations and against the possible 
uncertainties in our knowledge of the cosmological model. In 
a future work we will present a more detailed application to 
the ISW effect that takes into account more realistic physical 
conditions for the galaxy populations. Finally, we would like 
to point out that the formalism can be easily generalised to 
include information from several galaxy surveys at the same 
time. 
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